Charlottesville CDC Places
Lee LeBoeuf
10/7/2021
Data Source
- These data come from the Center for Disease Control and Prevention (CDC) PLACES: Local Data for Better Health.
- More information about the CDC PLACES data can be found here.
- Data shown here were downloaded using the CDC PLACES API in October of 2021. Documentation for the API can be found here and documentation for the data can be found here.
The Data Used Here
- PLACES is a collaboration between the CDC, the Robert Wood Johnson Foundation, and the CDC Foundation.
- The data provide model-based estimates for chronic disease risk factors, health outcomes, and clinical prevention services. According the the website, "CDC uses an innovative peer-reviewed multilevel regression and poststratification (MPR) approach that links geocoded health surveys and high spatial resolution population demographic and socioeconomic data."
- CDC pulls data from the CDC Behavioral Risk Factor Surveillance System, the Census 2010 population, and the American Community Survey estimates.
- The full data from the CDC include 27 measures: 5 unhealthy behaviors, 13 health outcomes, and 9 prevention practices. The data used here include 9 health outcomes, 3 prevention practices, and 5 unhealthy behaviors. All groupings shown below are consistent with the CDC groupings of each outcome.
- All of the data values represent the predicted percent prevalence of a given outcome within a census tract. Data for some outcomes are from 2017, but most are from 2018.
- Data are available at the tract levels for Albemarle, Charlottesville, Louisa, Greene, Fluvanna, Nelson counties.
Variable descriptions
These summaries are based on data at the block group level from the years 1980-2020. Data are also available at the block, tract, and county levels.
meta %>%
filter(su_tract == 1) %>%
select(varname, about) %>% as.list()## $varname
## [1] "locationname" "countyname"
## [3] "totalpopulation" "Coronary_Heart_Disease2018"
## [5] "Binge_Drinking2018" "Mental_Health2018"
## [7] "High_Blood_Pressure2017" "Physical_Inactivity2018"
## [9] "Diabetes2018" "Current_Smoking2018"
## [11] "Cancer_except_skin2018" "Current_Asthma2018"
## [13] "Dental_Visit2018" "High_Cholesterol2017"
## [15] "COPD2018" "Obesity2018"
## [17] "Physical_Health2018" "Health_Insurance2018"
## [19] "less_than_sevenhr_sleep2018" "Annual_Checkup2018"
##
## $about
## [1] "Tract FIPS code"
## [2] "County name"
## [3] "Total population of the tract as of the 2010 census"
## [4] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had angina or coronary heart disease. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [5] "Adjusted percent of respondents aged >= 18 years who report having five or more drinks (men) or four or more drinks (women) on an occassion in the past 30 days. Not an indicator of the frequency of binge drinking or the specific amount of alcohol consumed."
## [6] "Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their mental health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
## [7] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have high blood pressure. Women who were told high blood pressure only during pregnancy and those who wree told they had borderline hypertension were not included. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate. Not a measure of current high blood pressure."
## [8] "Adjusted percent of respondents aged >= 18 years who answered \"no\" to the following question: \"During the past month, other than your regular job, did you participate in any physical activity or exercise such as running, calisthenics, golf, gardening, or walking for exercise?\" Only captures information about nonoccupational physical activity."
## [9] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had diabetes (other than diabetes during pregnancy). Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [10] "Adjusted percent of respondents aged >= 18 years who report having smoked >= 100 cigarettes in their lifetime and currently smoke every day or some days. Not a measure of lifetime or current number of cigarettes smoked, and each of these factors can affect the risk for acquiring chronic disease from smoking. Not not measure intent or attempts to quit smoking among smokers or exposure to secondhand smoke among nonsmokers."
## [11] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have any type of cancer except skin cancer. Not specific to cancer type. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [12] "Adjusted percent of respondents aged >= 18 years who answer \"yes\" to both of the following questions: (1) \"Have you ever been told by a doctor, nurse, or other health professional that you have asthma?\" and (2) \"Do you still have asthma?\" This indicator requires doctor diagnosis, which may not include all persons with asthma"
## [13] "Adjusted percent of respondents aged >= 18 years who report having been to the dentist or dental clinic in the previous year."
## [14] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had high cholesterol out of all respondents who have had their cholesterol checked in the last 5 years. Reliability can be low on this indicator because patients might not know what specific tests have been performed on their blood samples unless their is an issue and many patients cannot afford to have their cholesterol checked."
## [15] "Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had COPD, emphysema, or chronic bronchitis. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
## [16] "Adjusted percent of respondents aged >= 18 years who have a BMI >= 30 kg/m^2 calculated from self-reported weight and height excluding respondents who were <3ft tall or >= 8ft; weighed <50lbs or >= 650 lbs; BMI < 12 or >= 100; pregnant women. Self-reports of height and weight lead to lower BMI estimates compared to height and weight measurements."
## [17] "Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their physical health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
## [18] "Adjusted percent of respondents aged 18-64 who report having no current health insurance coverage. All persons >= 65 are eligible for Medicare. This indicator is likely to be an underestimate due to how variable health insurance coverage can be."
## [19] "Adjusted percent of respondents aged >= 18 years who report usually getting insufficient sleep (<7 hours for those aged >= 18 years, on average, during a 24-hour period). Indicator does not measure variations in sleep duration, quality of sleep, or specific sleep problems."
## [20] "Adjusted percent of respondents aged >= 18 years who report having been to a doctor for a routine checkup (e.g., a general physical exam, not an exam for specifiic injury, illness, condition) in the previous year."glimpse(cvldat)## Rows: 50
## Columns: 20
## $ locationname <dbl> 51540000600, 51540000800, 51540001000, 515…
## $ countyname <chr> "Charlottesville", "Charlottesville", "Cha…
## $ totalpopulation <int> 4351, 3642, 2783, 3324, 3305, 2791, 4278, …
## $ Coronary_Heart_Disease2018 <dbl> 2.2, 7.1, 5.1, 5.8, 5.9, 3.5, 3.9, 3.9, 3.…
## $ Binge_Drinking2018 <dbl> 19.2, 15.8, 19.1, 16.7, 16.2, 21.2, 20.7, …
## $ Mental_Health2018 <dbl> 21.1, 11.2, 10.3, 16.2, 18.2, 15.0, 12.9, …
## $ High_Blood_Pressure2017 <dbl> 13.6, 33.7, 27.8, 33.5, 33.4, 21.6, 20.6, …
## $ Physical_Inactivity2018 <dbl> 23.1, 21.2, 16.1, 27.2, 29.8, 18.8, 16.2, …
## $ Diabetes2018 <dbl> 4.0, 10.5, 7.9, 12.1, 12.8, 6.7, 5.9, 6.6,…
## $ Current_Smoking2018 <dbl> 19.1, 14.2, 11.5, 21.8, 24.5, 15.4, 12.2, …
## $ Cancer_except_skin2018 <dbl> 1.4, 8.4, 6.9, 5.0, 4.4, 3.6, 4.7, 5.4, 1.…
## $ Current_Asthma2018 <dbl> 10.4, 8.7, 8.0, 10.7, 11.3, 9.5, 8.8, 8.1,…
## $ Dental_Visit2018 <dbl> 56.6, 73.2, 77.6, 60.0, 54.1, 68.9, 75.6, …
## $ High_Cholesterol2017 <dbl> 14.9, 36.1, 31.4, 31.0, 30.1, 23.3, 25.6, …
## $ COPD2018 <dbl> 4.0, 6.7, 4.5, 7.0, 7.6, 4.2, 4.1, 4.0, 5.…
## $ Obesity2018 <dbl> 24.1, 27.6, 26.2, 35.9, 37.7, 26.5, 23.2, …
## $ Physical_Health2018 <dbl> 9.6, 11.3, 8.7, 13.8, 15.2, 8.9, 8.0, 7.8,…
## $ Health_Insurance2018 <dbl> 19.4, 12.7, 9.8, 18.6, 21.1, 14.1, 10.8, 1…
## $ less_than_sevenhr_sleep2018 <dbl> 40.2, 34.4, 33.1, 42.3, 44.3, 36.9, 33.7, …
## $ Annual_Checkup2018 <dbl> 72.1, 80.5, 78.1, 78.2, 77.7, 74.7, 75.8, …cvldat %>% select(Coronary_Heart_Disease2018, Binge_Drinking2018, Mental_Health2018, High_Blood_Pressure2017, Physical_Inactivity2018, Diabetes2018, Current_Smoking2018,
Cancer_except_skin2018, Current_Asthma2018, Dental_Visit2018, High_Cholesterol2017, COPD2018, Obesity2018, Physical_Health2018, Health_Insurance2018, less_than_sevenhr_sleep2018,
Annual_Checkup2018) %>%
select(where(~is.numeric(.x))) %>%
as.data.frame() %>%
stargazer(., type = "text", title = "Summary Statistics", digits = 2,
summary.stat = c("mean", "sd", "min", "median", "max"))##
## Summary Statistics
## =============================================================
## Statistic Mean St. Dev. Min Median Max
## -------------------------------------------------------------
## Coronary_Heart_Disease2018 5.92 1.78 1.30 6.00 9.00
## Binge_Drinking2018 17.22 1.89 14.40 16.75 21.50
## Mental_Health2018 13.55 3.92 8.30 12.80 29.70
## High_Blood_Pressure2017 30.03 6.89 11.00 31.75 38.50
## Physical_Inactivity2018 21.54 4.59 14.00 21.25 31.60
## Diabetes2018 9.36 2.73 3 9.4 14
## Current_Smoking2018 16.52 4.37 8.70 15.25 24.50
## Cancer_except_skin2018 6.58 2.11 0.80 7.20 9.20
## Current_Asthma2018 9.28 1.14 7.40 9.15 13.20
## Dental_Visit2018 70.31 7.92 46.20 71.30 83.70
## High_Cholesterol2017 32.01 6.92 10.70 34.25 39.80
## COPD2018 6.43 1.80 3 6.2 11
## Obesity2018 30.14 5.16 20.10 28.40 40.60
## Physical_Health2018 11.49 2.42 7.10 11.35 16.80
## Health_Insurance2018 13.72 3.78 6.70 13.55 23.80
## less_than_sevenhr_sleep2018 35.75 3.45 28.70 35.65 44.30
## Annual_Checkup2018 78.33 2.67 71.50 79.05 82.10
## -------------------------------------------------------------Visual distribution
These distributions show data at the block group level from the years 1980-2020. Data are also available at the block, tract, and county levels.
longdat <- cvldat %>% select(c(locationname, Coronary_Heart_Disease2018, Binge_Drinking2018, Mental_Health2018, High_Blood_Pressure2017, Physical_Inactivity2018, Diabetes2018, Current_Smoking2018,
Cancer_except_skin2018, Current_Asthma2018, Dental_Visit2018, High_Cholesterol2017, COPD2018, Obesity2018, Physical_Health2018, Health_Insurance2018, less_than_sevenhr_sleep2018,
Annual_Checkup2018)) %>% pivot_longer(-locationname, names_to = "measure", values_to = "value")
longdat$measure <- factor(longdat$measure,
levels = c("Coronary_Heart_Disease2018", "Binge_Drinking2018", "Mental_Health2018", "High_Blood_Pressure2017", "Physical_Inactivity2018", "Diabetes2018", "Current_Smoking2018", "Cancer_except_skin2018", "Current_Asthma2018", "Dental_Visit2018", "High_Cholesterol2017", "COPD2018", "Obesity2018", "Physical_Health2018", "Health_Insurance2018", "less_than_sevenhr_sleep2018", "Annual_Checkup2018"))
longdat %>%
ggplot(aes(x = value, fill = measure)) +
scale_fill_viridis(option = "Blues", discrete = TRUE, guide = FALSE) +
geom_histogram() +
facet_wrap(~measure, scales = "free", ncol = 4)## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.## Warning in viridisLite::viridis(n, alpha, begin, end, direction, option): Option
## 'Blues' does not exist. Defaulting to 'viridis'.
meta %>%
filter(varname %in% c("Coronary_Heart_Disease2018", "Binge_Drinking2018", "Mental_Health2018", "High_Blood_Pressure2017", "Physical_Inactivity2018", "Diabetes2018", "Current_Smoking2018", "Cancer_except_skin2018", "Current_Asthma2018", "Dental_Visit2018", "High_Cholesterol2017", "COPD2018", "Obesity2018", "Physical_Health2018", "Health_Insurance2018", "less_than_sevenhr_sleep2018", "Annual_Checkup2018")) %>%
mutate(label = paste0(varname, ": ", about)) %>%
select(label) %>%
as.list()$label [1] "Coronary_Heart_Disease2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had angina or coronary heart disease. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[2] "Binge_Drinking2018: Adjusted percent of respondents aged >= 18 years who report having five or more drinks (men) or four or more drinks (women) on an occassion in the past 30 days. Not an indicator of the frequency of binge drinking or the specific amount of alcohol consumed."
[3] "Mental_Health2018: Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their mental health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
[4] "High_Blood_Pressure2017: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have high blood pressure. Women who were told high blood pressure only during pregnancy and those who wree told they had borderline hypertension were not included. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate. Not a measure of current high blood pressure." [5] "Physical_Inactivity2018: Adjusted percent of respondents aged >= 18 years who answered "no" to the following question: "During the past month, other than your regular job, did you participate in any physical activity or exercise such as running, calisthenics, golf, gardening, or walking for exercise?" Only captures information about nonoccupational physical activity."
[6] "Diabetes2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had diabetes (other than diabetes during pregnancy). Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[7] "Current_Smoking2018: Adjusted percent of respondents aged >= 18 years who report having smoked >= 100 cigarettes in their lifetime and currently smoke every day or some days. Not a measure of lifetime or current number of cigarettes smoked, and each of these factors can affect the risk for acquiring chronic disease from smoking. Not not measure intent or attempts to quit smoking among smokers or exposure to secondhand smoke among nonsmokers."
[8] "Cancer_except_skin2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they have any type of cancer except skin cancer. Not specific to cancer type. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[9] "Current_Asthma2018: Adjusted percent of respondents aged >= 18 years who answer "yes" to both of the following questions: (1) "Have you ever been told by a doctor, nurse, or other health professional that you have asthma?" and (2) "Do you still have asthma?" This indicator requires doctor diagnosis, which may not include all persons with asthma"
[10] "Dental_Visit2018: Adjusted percent of respondents aged >= 18 years who report having been to the dentist or dental clinic in the previous year."
[11] "High_Cholesterol2017: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had high cholesterol out of all respondents who have had their cholesterol checked in the last 5 years. Reliability can be low on this indicator because patients might not know what specific tests have been performed on their blood samples unless their is an issue and many patients cannot afford to have their cholesterol checked."
[12] "COPD2018: Adjusted percent of respondents aged >= 18 years who report ever having been told by a health professional that they had COPD, emphysema, or chronic bronchitis. Based on being diagnosed and respondent recall of diagnosis, so might be underestimate."
[13] "Obesity2018: Adjusted percent of respondents aged >= 18 years who have a BMI >= 30 kg/m^2 calculated from self-reported weight and height excluding respondents who were <3ft tall or >= 8ft; weighed <50lbs or >= 650 lbs; BMI < 12 or >= 100; pregnant women. Self-reports of height and weight lead to lower BMI estimates compared to height and weight measurements."
[14] "Physical_Health2018: Adjusted percent of respondents aged >= 18 years who report 14 or more days during the past 30 days during which their physical health was not good. Based self-assessment only and does not have an objective health component, so it's difficult to assess reliability and validity."
[15] "Health_Insurance2018: Adjusted percent of respondents aged 18-64 who report having no current health insurance coverage. All persons >= 65 are eligible for Medicare. This indicator is likely to be an underestimate due to how variable health insurance coverage can be."
[16] "less_than_sevenhr_sleep2018: Adjusted percent of respondents aged >= 18 years who report usually getting insufficient sleep (<7 hours for those aged >= 18 years, on average, during a 24-hour period). Indicator does not measure variations in sleep duration, quality of sleep, or specific sleep problems."
[17] "Annual_Checkup2018: Adjusted percent of respondents aged >= 18 years who report having been to a doctor for a routine checkup (e.g., a general physical exam, not an exam for specifiic injury, illness, condition) in the previous year."
Health outcomes
Asthma
pal <- colorNumeric("Blues", domain = mapdat$Current_Asthma2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Current_Asthma2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with asthma: ", mapdat$Current_Asthma2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Current_Asthma2018),
title = "Predicted % of adults <br> with current asthma <br> in 2018", opacity = 0.7)Coronary Heart Disease
pal <- colorNumeric("Blues", domain = mapdat$Coronary_Heart_Disease2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Coronary_Heart_Disease2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with Coronary Heart Disease: ", mapdat$Coronary_Heart_Disease2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Coronary_Heart_Disease2018),
title = "Predicted % of adults <br> with Coronary Heart Disease <br> in 2018", opacity = 0.7)Cancer (excluding skin cancer)
pal <- colorNumeric("Blues", domain = mapdat$Cancer_except_skin2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Cancer_except_skin2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with cancer: ", mapdat$Coronary_Heart_Disease2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Cancer_except_skin2018),
title = "Predicted % of adults <br> with cancer (excluding skin) <br> in 2018", opacity = 0.7)Diabetes
pal <- colorNumeric("Blues", domain = mapdat$Diabetes2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Diabetes2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults with diabetes: ", mapdat$Diabetes2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Diabetes2018),
title = "Predicted % of adults <br> with diabetes <br> in 2018", opacity = 0.7)Prevention
Lack of health insurance
pal <- colorNumeric("Blues", domain = mapdat$Health_Insurance2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Health_Insurance2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults without health insurance: ", mapdat$Health_Insurance2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Health_Insurance2018),
title = "Predicted % of adults <br> without health insurance <br> in 2018", opacity = 0.7)Routine check-up in the last year in 2018
pal <- colorNumeric("Blues", domain = mapdat$Annual_Checkup2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Annual_Checkup2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who went <br> to the doctor in last year: ", mapdat$Annual_Checkup2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Annual_Checkup2018),
title = "Predicted % of adults <br> who went to the doctor <br> in the last year <br> in 2018", opacity = 0.7)Visits to dentist/dental clinic in the last year in 2018
pal <- colorNumeric("Blues", domain = mapdat$Dental_Visit2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Dental_Visit2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who went <br> to the dentist in last year: ", mapdat$Dental_Visit2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Dental_Visit2018),
title = "Predicted % of adults <br> who went to the dentist <br> in the last year <br> in 2018", opacity = 0.7)Unhealthy Behaviors
Curent smoking in 2018
pal <- colorNumeric("Blues", domain = mapdat$Current_Smoking2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Current_Smoking2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who smoke: ", mapdat$Current_Smoking2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Current_Smoking2018),
title = "Predicted % of adults <br> who currently smoke <br> in 2018", opacity = 0.7)Obesity
pal <- colorNumeric("Blues", domain = mapdat$Obesity2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$Obesity2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who have a BMI >= 30: ", mapdat$Obesity2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$Obesity2018),
title = "Predicted % of adults <br> who had a BMI >= 30 <br> in 2018", opacity = 0.7)Less than 7 hours of sleep
pal <- colorNumeric("Blues", domain = mapdat$less_than_sevenhr_sleep2018)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data = mapdat,
fillColor = ~pal(mapdat$less_than_sevenhr_sleep2018),
weight = 1,
opacity = 1,
color = "white",
fillOpacity = 0.6,
highlight = highlightOptions(
weight = 2,
fillOpacity = 0.8,
bringToFront = T),
popup = paste0("GEOID: ", mapdat$GEOID, "<br>",
"Percent of adults who average <7 hrs <br> of sleep a night in 2018: ", mapdat$less_than_sevenhr_sleep2018)) %>%
leaflet::addLegend("bottomright", pal = pal, values = (mapdat$less_than_sevenhr_sleep2018),
title = "Predicted % of adults <br> who average <7 hrs <br> of sleep a night in 2018", opacity = 0.7)